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Abstract. In this work, a statistical analysis of the distribution of daily fluctuations of the IPC, the 
Mexican Stock Market Index is presented. A sample of the IPC covering the 13-year period 04/19/1990 - 
08/21/2003 was analyzed and the cumulative probability distribution of its daily logarithmic variations 
studied. Results showed that the cumulative distribution function for extreme variations, can be described 
by a Pareto-Levy model with shape parameters a = 3.634 ± 0.272 and a = 3.540 ± 0.278 for its positive 
and negative tails respectively. This result is consistent with previous studies, where it has been found 
that 2.5 < a < 4 for other financial markets worldwide. 
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Presentamos un analisis estadfstico de la distribution de fluctuaciones diarias del mdice de la Bolsa Mexi- 
cana de Valores, el llamado IPC (Indice de Precios y Cotizaciones). Estudiamos la funcion de distribution 
acumulativa de las diferencias logarftmicas diarias calculadas a partir de una muestra del IPC que cubre 
un periodo de 13 anos, que empieza el 19/04/1990 y finaliza el 21/08/2003. Hallamos que esta funcion de 
distribution acumulativa puede describirse para los valores extremos de estas diferencias mediante una 
distribution de Pareto-Levy (ley potencia) con exponentes a = 3.634 ± 0.272 y a — 3.540 ± 0.278 en sus 
colas positiva y negativa respectivamente. Este resultado es consistente con estudios previos que muestran 
que 2.5 < a < 4 para los mercados financieros de diferentes partes del mundo. 

PACS. 01.75.+m Science and Society - 02.50.-r Probability theory, stochastic processes and statistics - 
89.65.Gh Economics; econophysics, financial markets, business and management - 89.90.+n Other areas 
of general interest to physicists 

1 Introduction Econophysics and it is becoming an emergent branch of 

physics by itself [T1l2TI!fll4TlolIu] . 

The behavior of extreme variations of economics indexes, In order to describe the behavior of the distribution of 
stock prices or even currencies, has been a topic of in- financial time series variations, several models have been 



terest in finance and economics, and its study becomes proposed. Some of them are: 

' relevant in the context of risk management and Financial „ . ,. , .. w=, 

- ' t>- i rrn TT i-i i u j-£n u ~~ Gaussian distribution 171. 

Risk lneory. However, these analyses are usually difficult T ,. , ,r* , . „ . 

r , n i r i — Log-gaussian distribution (Geometric fJrowman Mo- 

to perform due to the small number of extreme obser- tion) El 

vations in the tails of the distributions of financial time „. ,', ^ ' v >_ -i ^ nrnrrn 

. . „ .i ,i ■ rn , • — Stable Levy distribution 9,10,11. 

series variations. Recently, the interest ol the physics com- „ , \ , „. , lo'i^mi 

• £ c -i i ± \ . i - Truncated Lew Distribution H2UT 314 . 

mum ty m the behavior ol financial markets, has strongly „ . ... ,. , . rrpirrTT — — 

i i , , ,i -i i-v. r , — Foisson like distribution 1151161 . 

increased, boosted by the availability of worldwide, elec- „ , — , ,, ,. 1 ., , 1 „ n-m , . 

. . ' . . ~ • • • i-rr — Power law Distribution with a ~ 3 1171 (Asymptoti- 



cally). 



tronically recorded financial data, giving rise to a differ- 
ent approach to confront problems arising from economics. 

The collection of methods and techniques originally devel- In this paper; a statistica i ana l ysis f the distribution of 

oped m the area of physics, which are currently applied daily variations of the IPC i is presented. It is organized as 

in the study of financial complex systems, is now called follows . In the remaining of this section, we briefly review 

a e-mail: hcoronel@uv.mx 1 The Mexican Stock Market index or Indice de Precios y 

b e-mail: alhernandez@uv.mx Cotizaciones for its Spanish meaning. 
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the Pareto-Levy distribution, and some of the phenom- 
ena it describes are mentioned. In next section, a very 
short introduction to the variables of common use in the 
study of financial index and prices variations is given. In 
section 3, we introduce the sample data analyzed in this 
work and some important statistical properties of financial 
time series variations (fat tails, clustering volatility, etc) 2 
are discussed, all the above in the context of the IPC data. 
In section 4 we explain and justify the procedure to esti- 
mate the Pareto-Levy exponent from data and we show 
the results concerning the fit on the tails of the cumula- 
tive distribution of the IPC daily logarithmic variations. 
Finally, last section is devoted to compare our results with 
some other related studies of different international stock 
markets previously reported. 



1.1 Pareto-Levy Distribution. Stable Distributions 

At this point, it is convenient to make a review of the 
definition of the Pareto-Levy distribution: 
An absolutely continuous random variable Y, is said to 
follow a Pareto-Levy distribution with parameters a and 
7, if its cumulative distribution function F has the form: 



F(y) := P{Y < y t } = 1 - 



= 1 - 



J_ 



(i) 



with yi > j/o j Vo — 7 an d oi > 0. When the condition 
a > 2 holds, the mean and variance of Y are both finite 
and by the central limit theorem, the sum of indepen- 
dent Pareto-Levy distributed random variables, converges 
in probability to a gaussian law. On the other hand, when 
a < 2, the Pareto-Levy distribution has infinite variance, 
and it is said that the distribution is stable. For a math- 
ematical treatment of these topics, consult 19 and |2()j . 
For a review from an econophysics point of view, refer- 
ences |oTl21| are recommended. 

The Pareto-Levy distribution is often known as the Power- 
law distribution, and its role in Physics and other ar- 
eas seems ubiquitous. In particular, we can illustrate this 
point by the following examples taken from finance: 

— N^t ) the number of trades in a given interval of time 
At, follows a power-law distribution with exponent 

— Pareto-Levy tails, with a ~ 3 for extreme variations 
of individual stocks prices [231 an d also for indexes of 
different leading stocks markets |17U24U25j . 

— Decay of volatility correlations follows a power law dis- 
tribution |25IIT7] . 

— The tail behavior of the cumulative distribution func- 
tion of volatility, is consistent with a power law distri- 
bution with exponent ~ 3 |28|. 



2 In finance, to these statistical properties of financial time 
series, jointly with other such as the properties of intermit- 
tency, asymmetry in time scales, absence of autocorrelations, 
gain/loss asymmetry etc., are called "Stylized facts" |18) . 



All the above suggests universality in financial complex 
systems and in order to explain the above facts, new mod- 
els and even theories, are currently being proposed 26 , 29 

LUSH. 



2 Study of variations of financial time series 

In the study of price variations of financial assets, many 
observables can be analyzed 0. If Y(t) is the value of the 
index at time t, some commonly used observables are: 

— Prices or indexes changes themselves, for some interval 
of time At : 



Z(t) :=Y{t + At) -Y(t) 
Deflated prices and index changes: 

Z D (t) := Z(t) x D(t) 



(2) 



(3) 



Where D{t) is a statistical factor or index called a discount 
or a deflation factor, and used to adjust the time value of 
money, enabling the comparison of prices while accounting 
for inflation, devaluation, etc. in different time periods. 

— Returns, defined as: 



R{t) 



Y(t + At) - Y(t) 



(4) 



— Differences of the natural logarithm of prices 3 , defined 
for some interval of time At as: 

S(t) := In Y(t + At) - In Y(i) (5) 

Each one having its own merits and disadvantages [H]. 
In this analysis, we have used the former variable S(t). 



3 Data sample and IPC variations 

The database containing the IPC series analyzed here is 
available at [221 and covers the 13-year period 04/19/1990 
- 08/21/2003. Figure U] shows the IPC evolution for this 
time period. We have used in our analysis the daily closure 
values of the IPC, that is, its recorded value at the end of 
each trading day 

In this work, our observable is S(t) as defined in equa- 
tion and we studied the tail behavior of P(S(t)) = 
1 - F{S(t)), for t = l,...,iV, where N = 3337 is our 
sample size and At = 1 day. 

Figure [2^ shows the histogram of Z(t), the IPC daily 
changes. It is interesting that this strongly symmetric and 
leptokurtic (fat tailed) distribution does not follow any 
well known model, which could appropriately describe the 
probability of events in its central region and in its tails 
at once. 



3 Note: Some authors call returns to the difference of the nat- 
ural logarithm of prices, and call normalized returns to our re- 
turns. For the case of high frequency data, each one approaches 
the other. 
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Fig. 1. IPC development for the 13 year period 04/19/1990- 
08/21/2003. 



Daily Diff Distribution 




Fig. 2. IPC variations, a) Histogram of Z(t), the daily changes 
of the IPC index. The distribution shows symmetry and fat 
tails, b) Histogram of S(t). Regions studied in this paper, 
S(t) > 0.035 and S(t) < -0.035, are indicated by two ver- 
tical lines, c) Same as above, but with a vertical logarithmic 
scale. The broken line corresponds to a gaussian density with 
the same mean and standard deviation than the S(t) distribu- 
tion. 



volatility 4 alternated with phases of a relative financial 
calm. This is a characteristic of financial time series called 
clustered volatility. 

The clustering phenomenon has not been completely un- 
derstood, however some successful models, such as the 
GARCH and Stochastic Volatility models have been pro- 
posed 5 PIP], 

In order to appreciate the magnitude of these strong fluc- 
tuations, gaussian values with the same mean and stan- 
dard deviation than those of the S(t) distribution were 
simulated. In the simulation, nearly no clustering is ap- 
preciated, as is shown in figure Oh- 
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Fig. 3. a) S(t) behavior for the period of time under study. 
Large variations in S(t), some of them as far as eight standard 
deviations from its mean, can be appreciated. It can also be 
seen that large variations tend to form clusters in time; this 
phenomena is called clustered volatility, b) Gaussian simula- 
tion already shown as a broken line in figure |5J;. Clustering is 
virtually not present. 



4 Parameter estimation from IPC empirical 
data 



Figures [2b and|2t show the distribution of our observable 
S(t). Figure^ shows that the distribution of S(t) appears 
to follow a gaussian. This is discarded after observing fig- 
ure the same distribution with a vertical logarithmic 
scale. There are too many extreme events visible almost 
as far as ten standard deviations from its mean. To easily 
comparing with, figure^ broken line is a gaussian scaled 
to the amplitude of S(t) and with the same mean and 
standard deviation than those of the S(t) series. 

Figure [3 shows how IPC variations are distributed; how- 
ever it does not give any insight about the dynamics of the 
stochastic process that governs them. Evolution of the IPC 
daily log-differences is displayed in figure We observe 
that large variations are not uniformly distributed over 
time, and that is possible to distinguish phases of higher 



The procedure to estimate a, the Pareto-Levy exponent 
from empirical data is straightforward: Pareto-Levy tails 
behave as straight lines in a logarithmic plot. Then, af- 
ter performing a linear fit of log P(S(t)) on log S(t), the 
obtained slope gives us an estimate of the exponent a 
of the Pareto-Levy Distribution. Clearly a distribution 
whose tails do not behave linearly in a log-log plot, can 
not be properly described by the Pareto-Levy model. 



4 In finance, volatility is a relative measure of price move- 
ment during a given time period. It can be modeled by the 
standard deviation of stock price changes. Econophysicists usu- 
ally use absolute returns to model the dynamics of volatility. 

5 The 2003 Nobel Price of Economy was awarded jointly to 
Robert F. Engle for a related topic:" ...for methods of analyzing 
economic time series with time- varying volatility (ARCH)". 
Engle owns a M. S. in Physics by Cornell University (1966). 
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The fit was carried out using data available in the tails 
regions \S(t)\ > 0.035. Both tails are marked with vertical 
lines in figure [2J). Those regions were chosen simply by 
examining the set of points for which the corresponding 
log-log plot behaves linearly. P(S(t)) was then reduced to 
79 and 93 events in its negative and positive tails respec- 
tively. Note that, as is showed in figure 01 and in order to 
deal with the undefined logarithmic scale for the left tail 
of P(S(i)), we have used —S(t) in our analysis. 
It was found, that for these regions, the tails of P(S(i)) 
decay following a power law model. A straight line pro- 
vides a good fit for them in a logarithmic plot. Both fits 
are shown in figure 0] 
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Fig. 4. Linear fitted tails in a log-log plot of the cumulative 
distribution function P(S(t)) on S(t). Right image positive tail. 
Left image negative tail. Fitted parameters are shown. 

Table^shows the estimated parameters and the 95% con- 
fidence intervals obtained from the lineal regression fit for 
the negative and positive tails of P(S(t)). 



Fitted region 



S(t) > 0.035 
S(t) < -0.035 



3.634 ± 0.272 
3.540 ± 0.278 



Table 1. Fitted Parameters plus minus twice the standard 
error of the estimates, for positive and negative tails. 



Table 2. Pareto exponent for some international Stock Mar- 
kets. Daily data (d). High Frequency data (*). 



Market 


a (Right tail) 


a (Left tail) 


Data Period 


AMEX 


2.84 ±0.12 


2.73 ±0.14 


01/94-12/95* 


NASDAQ 








NYSE 








(Combined) 








(USA) 








El 








S& P 500 


3.66 ±0.011 


3.61 ±0.11 


1962-1996 d 


(USA) 


3.39 ± 0.05 


3.37 ±0.07 


1984-1996* 


EH 








DAX 


2.4 (minutely) 




10/97-12/99* 


(Germany) 


to 3.5 (hourly) 




1959-2001 d 










NIKKEI 


3.05 ±0.16 




1984-1997 d 


(Japan) 

EH 








Hang-Seng 


3.03 ±0.16 




1986-1997 d 


(Hong Kong) 
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of stock markets of emergent economies, it seems they may 
belong to a different universality class, some studies |3bl 
137) show that the return distributions from emergent mar- 
kets have fatter tails than the observed in developed mar- 
kets. 

In summary, it has been shown that the cumulative prob- 
ability of daily extreme logarithmic changes of the Mexi- 
can IPC index, can be approximated by the Pareto-Levy 
model, with exponents a = 3.634 ±0.272 and a = 3.540 ± 
0.278 for its positive and negative tails respectively. As 
a consequence of these values, we can affirm that the 
stochastic process that governs the time series S(t), is well 
outside the Levy stable regime (0 < a < 2). 

Acknowledgments 
A.R.H.M thanks Conacyt-Mexico for financial support pro- 
vided under Grant No. 44598F. Also wishes to thank to N. 
Bagatella, N. Cruz, A. Guerra, G. Mandujano and R.de la 
Mora for their useful suggestions. We specially appreciate 
the bibliographical support given by P. Giubellino. 
Analysis in this work was performed using ROOT |38j . 



5 Discussion 

Results shown in tableware consistent with similar stud- 
ies, where the Pareto-Levy model with 2.5 < a < 4.0, has 
been found be useful for describe the behavior of extreme 
variations of diverse financial markets. Table F2] summa- 
rizes results of some of these studies. 
High frequency studies of price variations, most of them 
performed for stock markets belonging to developed coun- 
tries show that the distribution of returns follows a Pareto- 
Levy form with exponent converging to a ~ 3 as At de- 
creases to time intervals of about one minute. For the case 
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